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ABSTRACT 

We present new results on the recently discovered 69 ms X-ray pulsar AXS J161730-505505, the sixth 
youngest sample of all known pulsars. We have undertaken a comprehensive X-ray observing campaign of 
AXS J161730-505505 with the ASCA, SAX, and XTE observatories and follow its long term spin-down history 
between 1989 and 1999, using these, archival GINGA and ASCA data sets, and the radio ephemeris. The spin- 
down is not simply described by a linear function as originally thought, but instead we find evidence of a giant 
glitch (| AP/P\ > 10~ 6 ) between 1993 August and 1997 September, perhaps the largest yet observed from a young 
pulsar. The glitch is well described by steps in P and P accompanied by a persistent P similar to those in the Vela 
pulsar. The pulse profile of AXS J 16 1730-505505 presents a single asymmetric peak which is maintained over all 
observation epochs. The energy spectrum is also steady over time, characterized by a highly absorbed power-law 
with a photon index T = 1 .4 ± 0.2, consistent with that found for other young rotation powered pulsars. 

Subject headings: pulsars: general — pulsars: individual (AXS J161730-505505, PSR J1617-5055) — X-rays: 
general — supernova remnant — stars: neutron 
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this object in archival soft X-ray images of this field suggested 
that the source might be highly absorbed, making it difficult to 
determine whether it is associated with a SNR. 



1. INTRODUCTION 

Radio pulsars are thought to be highly magnetized (~ 10 12 
G), rapidly spinning neutron stars whose luminosity is pow- 
ered by rotational energy loss. The study of young (< 10 5 yrs) 
rotation-powered pulsars provides an important laboratory for 
understanding the early evolution (thermal, spin, and magnetic) 
of these embers of stellar collapse. The most energetic of these 
pulsars are observable at X-ray wavelengths, which allow us to 
probe these extreme, but rare, examples. 

In addition to uniform spin-down corresponding to the ro- 
tational energy loss, rotation powered pulsars show sudden 
discontinuities in their rotation periods (see Lyne & Graham- 
Smith 1998, chapter 6). These rare phenomena, known as 
"glitches" are considered to arise from sudden changes in the 
configuration of super-dense material in the neutron star inte- 
rior. To date, a total of 71 glitches with \AP/P\ > 10~ 9 have 
been reported in 30 pulsars (Urama & Okeke 1999). Observa- 
tion of pulsar glitches gives us insights into the structure and 
physical processes inside the neutron stars, such as the interac- 
tions of neutron superfluid and crust components (e.g., Ander- 
son & Itoh 1975). 

In this article we report the detection of a giant glitch 
from AXS J 16 1730-505505 by using the newly acquired multi- 
mission X-ray data. This source is an unusual case of a young 
pulsar discovered first by its X-ray emission, revealed during 
the course of an archival X-ray study of the SNR RCW 103 
(Gotthel, Petre, & Hwang 1997). Further analysis detected 
highly significant pulsations from photons attributed to this 
source (Torii et al. 1998). The 69 ms pulse period was consis- 
tent with that reported from a GINGA observation of the region 
(Aoki, Dotani, & Mitsuda 1992). Recent radio observation has 
confirmed AXS J161730-505505 as a young energetic rotation- 
powered pulsar (Kaspi et al. 1998). The lack of evidence for 
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2. OBSERVATIONS 

A set of day-long X-ray observations of the field containing 
RCW 103, IE 161348-5055, and AXS J16 1730-505505 were 
performed with the ASCA (Tanaka, Inoue, & Holt 1994), SAX 
(Boellaefa/. 1997), and XTE (Bradt et al. 1993) observatories. 
We summarize in table 1 the set of observations presented in 
this work. 

Imaging data were acquired with the Gas Imaging Spec- 
trometers (GIS) on-board ASCA and with the Medium Energy 
Concentrators (MECS) instruments on SAX. These instruments 
have moderate imaging (~ 2') and spectral resolution (~ 8% 
at 6 keV) over an energy band pass of 0.7 - 10 keV (GIS) and 
1.5- 12 keV (MECS), with a field-of-view (FOV) large enough 
to cover the SNR and pulsar simultaneously. Non-imaging data 
were obtained with the Proportional Counter Array (PCA) on- 
board XTE which provides broader energy band-pass (2 — 60 
keV) at lower spectral resolution (~ 16% at 6 keV). The PCA 
FOV is roughly circular with a ~ 1° FWHM response. The 
GIS data were collected in the highest time resolution mode 
(0.5 ms or 61 /is, depending on data acquisition mode) whose 
measured absolute accuracy is 200 /is in this mode (Saito et al. 
1997). The PCA data were collected using the Good Xenon 
mode with 0.9/is timing resolution. For the current analysis, 
the absolute timing uncertainty is ~ 100/is (Rots et al. 1998). 
Photons collected by the MECS are time tagged with 15 /is res- 
olution. We do not include data from ASCA's Solid-state Imag- 
ing Spectrometers (SIS), as the pulsar fell just off the edge of its 
FOV. Nor do we include data from SAX's other instruments as 
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the observing time is insufficient for these instruments to mea- 
sure the pulsar periodicity reliably, as a part of the observation 
was interrupted prematurely. 

Each data set was processed through its standard pipeline re- 
duction for that mission and edited to exclude times of high 
background contamination using the standard screening crite- 
ria. This rejects time intervals of South Atlantic Anomaly pas- 
sages, Earth occultations, bright Earth limb in the FOV (ASCA 
and SAX only), and other periods of high particle activity. The 
resulting effective observation times are summarized in Table 
1. For each observation, event data from all detectors were 
co-added and the arrival times of each event were corrected 
to the solar system barycenter using the software TIMECONV 
(ASCA), BARYCONV (SAX), or FXBARY (XTE). 

3. RESULTS 
3.1. Timing 

The X-ray images obtained with both ASCA and SAX above 
3 keV reveal AXS J16 1730-505505 4' outside the SNR shell 
(see figure la and lb of Gotthelf et al. 1997; Gotthelf, Petre, & 
Vasisht 1999). To increase the signal-to-noise ratio for detect- 
ing pulsations from the pulsar, we extracted photons from an 
8' diameter aperture centered on the pulsar, restricting the en- 
ergy range of extracted photons to 3 - 1 ke V for GIS and 3-12 
keV for the other instruments. For the PCA data, in this energy 
band, we further restrict our search to Layer 1 data only, which 
provides the best sensitivity for a Crab-like spectrum; For the 
higher energy analysis afforded by the PCA, above ^12 keV, 
we used data from all three PCA layers. 

We searched each data set for the expected 69 ms period pre- 
dicted from the initial period and period derivative measure- 
ment (Torii et al. 1998). A periodgram was constructed using 
the x 2 statistic to test against a null hypothesis. For each trial 
period, we folded the data into 10 bins and computed the x 2 of 
the resultant profile. We search a narrow range of periods cen- 
tered on the expected period ±0.1 ms, sampled in increments 
of 0.1 x P 2 /T, where T is the observation duration, and P is the 
test period. A highly significant signal was detected from each 
of our data sets. 

As well as the newly obtained data, we have re-analyzed the 
previous GINGA and ASCA data (Aoki et al. 1992; Torii et al. 
1998) in a uniform way and revised the period and its error by 
using the method of Leahy (1987). Our X-ray timing results 
derived from these 13 measurements are listed in table 1, along 
with an updated radio ephemeris (Kaspi 1999, Private commu- 
nication). 

3.2. Spectrum 

We search for spectral dependence of the pulse profile by 
comparing the folded light curves in several energy bands. No 
strong energy dependence is evident in the energy resolved light 
curves. Furthermore, the pulse amplitude and pulse profile re- 
mained unchanged between observational epochs. 

We examined the ASCA and SAX data on 
AXS J 16 1730-505505 for any long term changes in its en- 
ergy spectrum or flux. As for the timing analysis, we re- 
strict our comparison to the energy range above 3 keV and 
extract photons from an 8' diameter aperture centered on the 
source. We fitted the spectrum with a power law function mod- 
ified by interstellar absorption. The absorption was fixed at 
6.8 x 10 22 cm -2 (Torii et al. 1998). Spectra from each ob- 
servation were found to be consistent with each other. Com- 



bining the 7 ASCA observations, we obtain the pulse phase 
averaged photon index, T = 1 .4 ± 0.2, and the observed flux of 
(3.6 ±0.2) x 10" 12 ergss _1 cm" 2 (90% confidence errors) in the 
3-10 keV range, which is consistent with the previous measure- 
ment (Torii et al. 1998). 

4. DISCUSSION 

A x 2 fit to the ai l 14 data points as summarized in table 1 
gives the mean spin-down rate of P = 1 .361 1(1) x 10" 13 s/s and 
P = 0.069347150(1) at MJD 50,000.0 (Figure 1). However, the 
quality of the fit is bad, with X 2 /d.o.f. = 4808/12. The residual 
of the fit shows a jump of AP ~-1.2x 10" 7 s between the ob- 
servations of 1993 August (MJD 49,217.6) and 1997 September 
(MJD 50,696.0) (Figure 2, top panel). Within the observation 
span of 10 years, the residual is neither periodic nor smooth. A 
sudden change in the period between MJD 49,217.6 and MJD 
50,696.0 is suggested. We consider the most likely explanation 
for these residuals is due to a glitch, similar to those observed 
in several young rotation-powered pulsars. 

Given the clear evidence of glitch activity, we next attempted 
to model the spin-down data with a single glitch followed by 
an exponential recovery. The data coverage is limited and we 
simply assumed the following relation for the spin-down. 

P(t) = P +P-(t-t ) + AP-exp(--^) (1) 

T 

where AP = for t < to and AP is a negative constant for 
t > to. This model contains five parameters, which, except 
for the depth of the glitch AP, are found to be independent 
of the time of the glitch, to. The derived parameters, Pq, P, 
AP, and r are summarized in table 2 for assumed values of 
t = 49, 300.0 MJD, t = 50, 000.0 MJD, and t = 50, 600.0 MJD. 
The residual for to = 50,000 is shown in the middle panel 
of figure 2. The quality of the fit is now characterized by 
X 2 /d.o.f. = 22.6/10. 

The size of the glitch depends strongly upon the unknown 
glitch epoch to- The fractional increase in rotation was found 
to be AP/P = -4.2 x 10" 6 for t = 50,000.0 MJD (fixed) while 
it changes between AP/P = -1 1 x 10" 6 for t = 49, 300.0 MJD 
(fixed) and AP/P = -1 .8 x 10" 6 for t = 50,600.0 MJD (fixed). 
The minimum fractional increase in rotation rate is therefore 
comparable with those of the largest known pulsar glitches 
(Lyne, et al. 1996a; Shemar & Lyne 1996). 

Using the above model, the recovery time following the 
glitch episode is found to be t = 700 days. This duration is 
somewhat unusual compared to the radio pulsars whose re- 
covery time is seen to bifurcate between r ~ 100 days and 
t > 1,000 days (Shemar & Lyne 1996). Our derived value 
lies squarely between these two time scales, perhaps due to the 
simple model we invoked which allows for only a single glitch. 
Because of the sparse data coverage between August 1993 and 
September 1997, however, we cannot determine if the recovery 
time could be expressed as the sum of the two timescales. 

Glitches found in radio pulsars may be classified into three 
groups (Lyne & Graham-Smith 1998). The first is a Crab-like 
glitch which is characterized by steps mainly in P. The second 
is a Vela-like glitch which is characterized by large changes 
in P (AP/P ~ -10~ 6 ) and the exponential recoveries. For the 
Vela pulsar, linear sawtooth changes in P have been observed 
between glitches (Lyne et al. 1996b). The third kind is often 
found in old pulsars which is characterized by a change in P. 

Since the large glitch found for AXS J 16 1730-505505 is 
similar in its magnitude to those in the Vela pulsar, a phe- 
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nomenological model taking into account the sawtooth behav- 
ior may be a good description. Apart from short term ef- 
fects, the spin-down for the Vela pulsar is expressed by a linear 
change in spin-down rate. Therefore, the following function is 
appropriate if the transient effects have already ceased by MJD 
50,696. 

P(t) = P + AP + (P + AP) ■ (t - 1 ) + l -P ■ (t - f ) 2 (2) 

Here, APrj = and AP = for t < to and they are constant val- 
ues (APo < and AP > 0) for t > to- This model contains 
six parameters. Again, the time of the glitch, to, had to be as- 
sumed. The derived parameters are summarized in table 3. For 
the condition that AP < 0, to was restricted to f ^ 50,205. 
The residual for fo = 50,000 is shown in the bottom panel of 
figure 2. The quality of the fit is significantly improved to 
X 2 /d.o.f. = 17.1/9. Compared to the fit to Equation (1), the 
f-test gives a chance probability of ~ 0.1. This result suggests 
that transient effects had indeed ceased by MJD 50,696. The 
reduced \ 2 is still larger than unity, suggesting the presence of 
timing noise and smaller glitches. 

We can estimate the expectancy of large (|AP/P| > 10" 7 ) 
glitches by using the semi-empirical relation based upon the 
superfiuid vortex unpinning model (Alpar & Baykal 1994). Us- 
ing this relation, the expected number of large glitches between 
the first GINGA observation and the last ASCA observations 
is 3.7. Therefore, there should have been about 4 glitches of 
|AP/P| > 10~ 7 . Indeed the residuals to the fit of Equation 
(1) or (2) still hint a small jump of AP = -5 x 10" 9 s be- 
tween MJD 51,263.8 (1999 Mar.) and MJD 51,394.2 (1999 



Aug). This may be another glitch of AP/P = -7 x 10~ 8 , much 
smaller in magnitude than the one near MJD 50,000, but still 
relatively large compared to those seen for most radio pulsars 
(Shemar & Lyne 1996). We have thus found a giant glitch of 
| AP/P | > 10" 6 at f ^ 50,000 MJD (49,218 < t < 50,696) 
and possibly a glitch of | AP/P| ~ 7 x 10" 8 at t ~ 5 1 , 300 MJD 
(51,264 <f < 51,394). 

The detection of a giant glitch from AXS J161730-505505 
gives a rare sample for studying the origin of pulsar glitches 
and the interior structure of neutron stars. In this context, reg- 
ular timing observation of this pulsar in the radio band is quite 
important to monitor the onset of a glitch and following tran- 
sient effects. Long term timing observations for measuring the 
braking index are desired. Searches for the highly absorbed X- 
ray emission from the supernova remnant associated with this 
pulsar should be undertaken with the next generation X-ray ob- 
servatories. 
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Table 1 

Timing Data for AXS J16 1730-505505. 



Number 


Satellite/ 


Instrument 


Obs Date" 


Exposure 


Epoch" 


Period" 




Observatory 






(s) 


(MJD) 


(ms) 


1 


GINGA 


LAC 


1989 Mar 4 


14,400 


47590.2 


69.31890(3) 


2 


ASCA 


GIS 


1993 Aug 17 


35,353 


49217.277625 


69.338019(2) 


3 


ASCA 


GIS 


1997 Sep 4 


58,063 


50696.736757 


69.355301(1) 


4 


XTE 


PCA 


1998 Jan 2 


33,712 


50816.132391 


69.356717(1) 


5 


Parkes 


64 m 


1998 Jan 15 




50829.7 


69.356889(7) d 


6 


SAX 


MECS 


1998 Sep 17 


16,288 


51073.892337 


69.359778(1) 


7 


ASCA 


GIS 


1999 Feb 12 


15,405 


51222.279184 


69.361532(4) 


8 


ASCA 


GIS 


1999 Feb 20 


16,232 


51229.708310 


69.361624(7) 


9 


XTE 


PCA 


1999 Mar 5 


37,253 


51243.247065 


69.3617824(5) 


10 


ASCA 


GIS 


1999 Mar 6 


20,700 


51243.354377 


69.361782(8) 


11 


SAX 


MECS 


1999 Mar 23 


38,197 


51261.461543 


69.362000(3) 


12 


ASCA 


GIS 


1999 Mar 26 


18,914 


51263.519288 


69.362021(3) 


13 


SAX 


MECS 


1999 Aug 4 


53,115 


51394.834928 


69.363568(2) 


14 


ASCA 


GIS 


1999 Sep 13 


17,796 


51434.716949 


69.364037(6) 



"Start time of the observation. 

"Middle time of the observation except for the radio ephemeris. 
c l<7 error of the last significant figure is shown in parentheses. 

d The radio ephemeris reported here (Kaspi 1999, Private communication) has been revised from that given 
inKaspiefa/. (1998). 



Table 2 

Best- fit Parameters of a Giant Glitch in AXS J161730-505505 for 

Equation (1). 



Parameter 




Values 




to [MJD] 


49,300 (fixed) 


50,000 (fixed) 


50,600 (fixed) 


Po [s]« 


0.069338991(2) 


0.069347220(3) 


0.069354272(5) 


P [s/s] 


1.3605(5) x 10~ 13 


1.3605(5) x 10~ 13 


1.3605(5) x 10~ 13 


AP [s] 


-8(2) x 10~ 7 


-2.9(2) x 10~ 7 


-1.26(4) x 10~ 7 


t [days] 


7(1) x 10 2 


7(1) x 10 2 


7(1) x 10 2 



"The period just before the glitch at to. 
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Table 3 

Best- fit Parameters of a Giant Glitch in AXS J161730-505505 for 

Equation (2). 



Parameter 




Values 




fo [MJD] 


49,300 (fixed) 


50,000 (fixed) 


50,205 (fixed) 


Pols] 


0.069338980(3) 


0.06934707(3) 


0.06934943(3) 


APo [s] 


-4.6(3) x 10~ 7 


-1.0(2) x 10~ 8 


-3.1(3) x 10~ 10 


P [s/s] 


1.344(3) x 10~ 13 


1.332(5) x 10~ 13 


1.328(5) x 10~ 13 


AP [s/s] 


5.8(7) x 10~ 15 


5.8(7) x 10~ 15 


5.8(7) x 10~ 15 


Pis-'] 


-2.0(3) x 10~ 23 


-2.0(3) x 10~ 23 


-2.0(3) x 10" 23 



CO 




CO 

CD L i I i i i i I i i i i I i i i i I i i i i I i i i i I i l 

° 4.7xl0 4 4.8xl0 4 4.9xl0 4 5xl0 4 5.1xi0 4 5.2xl0 4 

Time [MJD] 

Fig. 1 .— The pulse period evolution of AXS J 16 1730-505505 over 10 years. 
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Fig. 2. — Top Panel) Pulse period residual as fitted by a linear function. Middle Panel) Pulse period residual as fitted to Equation (1) 
with a glitch at MJD 50,000. Bottom Panel) Pulse period residual as fitted to Equation (2) with a glitch at MJD 50,000. 



